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ABSTRACT 

We present period-luminosity relations for more than 3,200 red variable stars in the Small 
Magellanic Cloud observed in the second phase of the Optical Gravitational Lensing Ex- 
periment (OGLE-II). Periods of multiply-periodic light curve solutions combined with the 
single-epoch 2MASS JHK$ magnitudes, reveal very similar distributions to those for the 
Large Magellanic Cloud in Paper I. The main features include four pulsating Asymptotic Gi- 
ant Branch (AGB) ridges, three distinct short-period sequences below the tip of the Red Giant 
Branch and two long-period sequences of ambiguous origin. We derive a relative distance 
modulus of the Clouds from the period-luminosity distributions for all stars of A/Zo = 0.44 
mag, which is in good agreement with recent independent results. The tip of the Red Giant 
Branch shows a colour and metallicity dependence that is in excellent agreement with the 
empirical results for globular clusters. We conclude that most variable stars below the TRGB 
are indeed RGB stars. 
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, ; 1 INTRODUCTION 



The period-luminosity (P-L) relations of long period variables on 
the AGB are well documented for the Large Magellanic Cloud. Ob- 
servations have revealed not only that Mira variables follow a tight 
near-infrared P-L relation (Glass & Lloyd Evans 1981, 2003; Feast 
et al. 1989), but also that semiregular variables pulsating in over- 
tone modes follow distinct and well-defined P-L relations (Wood et 
al. 1999, Wood 2000). In recent years, there has been strong interest 
in pulsating red giants, which is reflected in the increasing number 
of independent analyses of large photometric databases (Cioni et 
al. 2001, 2003, Noda et al. 2002, Lebzelter et al. 2002, Glass & 
Schultheis 2003, Kiss & Bedding 2003, Ita et al. 2003, Groenewe- 
gen 2004). 

An interesting by-product was the recognition of a distinct and 
very numerous group of red variables below the tip of the Red Gi- 
ant Branch (TRGB). First, Ita et al. (2002) noted that the luminosity 
function (LF) of red giant variable stars in the LMC shows a sharp 
feature at the expected brightness of the TRGB. Earlier, Alves et 
al. (1998) and Wood (2000) suggested that thermally pulsing AGB 
stars can mimick a two-peaked luminosity function, and thus the 
LF shape alone is not sufficient to establish the existence of TRGB 
variables. However, we have shown in Paper I (Kiss & Bedding 
2003) that the LMC period-luminosity relations for stars above and 
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below the TRGB show a relative period shift that is consistent with 
evolutionary temperature difference between RGB and AGB stars. 
This result has since been confirmed by Ita et al. (2003). The pos- 
sibility of observable pulsations in RGB stars therefore seems to 
be well established, providing a new area for asteroseismological 
considerations. In this Letter we provide further evidence for ex- 
tragalactic RGB variables through a comparative analysis of red 
variables in the Small Magellanic Cloud (SMC). 

Red giant pulsators in the SMC have a far less extensive liter- 
ature than those of the LMC. The largest pre-microlensing datasets 
for the SMC were published by Wood et al. (1981), Lloyd Evans et 
al. (1988) and Sebo & Wood (1994), all consisting of a few dozen 
variables at most. Microlensing surveys (MACHO, EROS, OGLE) 
have changed the situation by discovering thousands of red vari- 
ables. Based on the extensive new datasets, it has became possi- 
ble to compare the Magellanic Clouds and the Galactic Bulge in 
search for metallicity effects on red giant pulsations and P-L rela- 
tions. Cioni et al. (2003) examined MACHO light curves of 458 
LPVs in the SMC that were detected by the Infrared Space Obser- 
vatory. They found very similar P-L sequences to those in the LMC 
and confirmed the relative overabundance of carbon-rich Mira stars 
(as recognized first by Lloyd Evans et al. 1988). The most extensive 
analysis to date is that of Ita et al. (2003), who determined the domi- 
nant period for 2,927 variables in the SMC using OGLE-II data and 
constructed P-L relations using single-epoch SIRIUS JHK survey 
data. Beside finding the same distributions as those in Paper I, Ita 
et al. (2003) derived zero-point differences of the P-L relations of 
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about 0. 1 mag. Some differences were also noted recently by Cioni 
et al. (2003) for the Magellanic Clouds and by Glass & Schultheis 
(2003) for the Galactic Bulge and the Magellanic Clouds. The latter 
investigators found both slightly different slopes and different ab- 
solute magnitude ranges for the overtone P-L sequences. These in- 
teresting results depend on the assumed distances to each of the dif- 
ferent P-L datasets. Here we concentrate on the multi-wavelength 
luminosity functions to sort out the distance issue, and examine the 
implications concerning the RGB pulsations. 



2 DATA ANALYSIS 

The basis of our analysis is the catalogue of OGLE-II (Udalski et 
al. 1997) data spanning four years from 1997 to 2000 (Zebrun et 
al. 2001). Of the total 7 square degrees of the OGLE Magellanic 
Cloud fields, about 2.5 square degrees were observed in the SMC. 
The observations were analysed by the OGLE team using a modi- 
fication of the Difference Image Analysis (Alard & Lupton, 1998, 
Wozniak 2000), which is the most advanced method available for 
detecting small variations in crowded fields (see, e.g., Bonanos & 
Stanek 2003). 

Typical light curves in the SMC consist of 280 points dis- 
tributed over ~ 1,100 days. Compared to the LMC, where individ- 
ual datasets covered slightly more than 1,200 days in 400 points, 
this means a reduced accuracy in period determination. We con- 
sider our periods reliable up to about 400-500 days, for which at 
least two cycles were observed. 

The data were analysed in the same manner as in Paper I. To 
summarise, we performed the following steps: 

(i) We cross-correlated the 15,038 OGLE-II variable stars in the 
SMC with the 2MASS All-Sky Point Source Catalog 1 ; within a 
search radius of 1" we found 10,361 stars with JHKs magnitudes. 
The tight constraint on coordinates ensured that we identified the 
infrared counterparts unambiguously. 

(ii) After excluding duplicates (< 10), we selected red giants ac- 
cording to the J — K s colour. The applied threshold was J — K s > 
0.9 mag. We have also excluded all stars with Ks > 14 mag be- 
cause of the strongly increasing errors of 2MASS magnitudes be- 
low that limit. The final sample consisted of 3,898 stars. 

(iii) We performed a period search on OGLE-II light curves by 
an iterative Fourier analysis (see Paper I for details), until a four- 
component harmonic fit was calculated. We kept only those fre- 
quencies larger than 8 x 10~ 4 d _1 (~ 1/T b s ) and with semi- 
amplitudes larger than 5 mmag. 

(iv) The resulting database contains 3,260 stars with 10,009 pe- 
riods (and amplitudes and phases), the mean /-band magnitude and 
2MASS JHK S single-epoch magnitudes. 

Our sample is slightly larger than that of Ita et al. (2003) be- 
cause of the different selection criteria. A more relevant difference 
is that we calculated a set of periods for every star, while Ita et al. 
(2003) determined only one period per star. We found multiple peri- 
odicity for the majority of stars (82% of the full sample), which can, 
therefore, be identified as common behaviour in red giants. There 
must be, of course, a fraction of spurious or physically unrelated 
periods, which need a more sophisticated method for identification 
and exclusion from future analyses. 
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Figure 1. P-L relations for red giants in the SMC (10,009 periods for 3,260 
stars). 

3 DISCUSSION 

We present the resulting P-L diagram in FigQ At first glance, the 
shape and extension of distributions are very similar to those we 
found for the LMC: (i) there is a noticeable drop of point density at 
K s ~ 12.70 mag, located close to the tip of the RGB (Cioni et al. 
2000); (ii) above the tip, there are four pulsating AGB sequences, 
which in the LMC were identified by Wood (2000) as stars pulsat- 
ing in the fundamental (F), first (10), second (20) and third (30) 
overtone modes (he labelled them as sequences C, B and A which, 
however, does not account for the fact that there are four ridges 
clearly distinguishable in the OGLE-II data); (iii) below the tip, 
there are three distinct sequences in the short-period range (P < 60 
days), of which R2 and R3 are continuations of 20 and 30; there 
might be a slight horizontal shift, similar to that found in the LMC 
(Paper I), but the ridges are less well defined, due to the smaller 
number of stars, so the possible period shift between 30 - R3 and 
between 20 - R2 is hardly measurable. The sequence Ri, just as 
for the LMC, lies somewhere between the extrapolated 10 and F 
ridges; (iv) for longer periods, there are two sequences, Li and L2, 
in the same relative position to the other ridges as in the LMC. 
Their interpretation is ambiguous; Wood (2000) suggested Li (=E 
in his paper) contains eclipsing binary variables, while L2 (=D) 
may consist of pulsating stars with unknown excitation mechanism 
or non-spherical stars with rotationally induced variatons (Olivier 
& Wood 2003, Wood 2004). 

Although multiperiodicity is a common feature, there is a 
remarkable difference above and below the tip of the RGB: for 
K s < 12.7 mag, 95%, while for K s > 12.7, only 62% of the 
light curves resulted in multiple periods. This, however, only re- 
flects the fact that amplitudes below the TRGB strongly decreases 
towards fainter magnitudes (see Sect. 3.3), thus the fixed 5 mmag 
amplitude threshold may be too high for most of the periods (if they 
exist). In addition, we stress there is not any noticeable change be- 
tween the resulting P-L distributions when plotting every star just 
once using the strongest period or using all periods with amplitudes 
exceeding a certain threshold (or even plotting every star just once 
with the second strongest period). 
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Table 1. Relative distance moduli and extinctions. Aft is the mean vertical 
distance between the P-L images; AAj refers to the differential extinction 
in the X band, while A/iq is the dereddened relative distance modulus. 
AE(B - V) = 0.085 (Westerlund 1997) and AA V = 0.26 mag were 
assumed. The uncertainty in Afi is about ±0.05 mag. 
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3.1 The relative distance modulus of the Clouds 

In the case of such complex structures as the P-L distributions of 
red variables in the Magellanic Clouds, it is difficult to measure 
the vertical distance between the P-L ridges directly. It has been 
a more common approach to assume a relative distance modulus 
and overplot the relations with the purpose of comparison, either 
in a simplified functional form (see Fig. 1 1 in Glass & Schultheis 
2003) or showing the full observational data (Fig. 10 in Ita et al. 
2003). Here we follow a different approach to show that, assuming 
a negligible mean zero-point difference for the set of relations, one 
can get consistent relative distance moduli and extinctions of the 
Clouds in all four photometric bands (IJHKs), which will be used 
later to characterize the dips of the luminosity functions. 

Instead of a visual examination of overplotted data, we cal- 
culated the cross-correlation functions of smoothed images of the 
P-L ridges (being the sharpest in Ks and the least-defined in 7). 
For this, we converted the scatter diagrams (Fig. Q and Fig. 1 in 
Paper I) to 500 x 500 pixel images as follows. The x-axis cov- 
ered logP between 1.2 and 3.1 (3.8 x 10~ 3 dex/pixel), while the 
y-axis covered 5 magnitudes, where the exact range depended on 
the studied band (e.g., from 9 mag to 14 mag in Ks). This way 
one pixel in vertical direction corresponds to 0.01 mag. Every pixel 
"intensity" value was set to the total number of stars in a 21 x 21 
pixel mask, centered on the actual pixel. This masking worked as a 
boxcar smoothing procedure, which produced clear images of the 
P-L ridges. The cross-correlation functions of the IJHKs P-L im- 
ages were calculated as functions of Afi, and their maxima showed 
the vertical shifts that give the best overlap of the images. In or- 
der to minimize the effects of spurious periods, we kept only stars 
with log P < 2.6 (i.e. P < 1 year); in other words, we used only 
360 x 500 pixel sub-images. 

The resulting shifts are presented in the second column of Ta- 
ble Q There is a strong wavelength dependence, which can be ex- 
plained by the differential reddening. The average E(B — V) for 
the LMC and the SMC are 0.15 mag and 0.065 mag, respectively 
(Westerlund 1997), which can be translated to A Ax differential 
extinction by adopting Rv =3.1 and the relative extinction coef- 
ficients (A x /A v ) in Schlegel et al. (1998). 

As can be seen in Table Q the dereddened relative distance 
moduli agree well within the error bars (~ ±0.05 mag in every 
band). Similar agreement exists for recently published results; for 
instance, Cioni et al. (2000) determined Afi — 0.44 mag. On the 
other hand, Ita et al. (2003) adopted Afi = 0.49 mag, while Glass 
& Schultheis (2003) used 0.50 mag. In the rest of the paper we 
adopt 0.44 mag with an estimated ±0.03 mag random error. 




Figure 2. Stellar magnitude distributions for OGLE-II variables in the SMC 
(thick red line) and in the LMC (thin black line), normalized by the survey 
area. 




Figure 3. The second derivatives of the smoothed luminosity functions. The 
arrows point to the maxima for the LMC (solid lines) and the SMC (dashed 
lines). 



3.2 The luminosity functions and the tip of the RGB 

In Fig. [2] we plot luminosity functions in all four bands of the 
OGLE-II+2MASS database: the thin solid line shows the LMC data 
after correction for the different survey areas (4.5 square degrees 
for the LMC vs. 2.5 square degrees for the SMC); the thick solid 
line denotes the SMC LF. 

Both Clouds have double-peaked LFs with similar shapes and 
wavelength dependences. The exact overlap of the faint ends of 
the SMC and area-normalized LMC distributions in JHKs shows 
that there is a severe sensitivity cut-off for J > 14.2 mag, H > 
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Table 2. Summary of TRGB magnitude values (the uncertainty is about 
0.04 mag). C2000 refers to TRGB values in Cioni et al. (2000). 
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13.5 and Ks > 13.3 mag. On the other hand, the difference for 
I suggests the drop of the OGLE-II detection efficiency occurs at 
fainter magnitudes, so that the faint limit in our sample is set by 
2MASS and not by OGLE. 

An interesting point is the changing shape of the AGB peak. It 
is the narrowest in / (where the full-width-at-half-maximum is 0.7 
mag) and going toward the redder bands, the stars spread over wider 
and wider luminosity ranges (the FWHM increases to 2 mag). This 
can be understood qualitatively in terms of the larger bolometric 
correction for the brightest and reddest stars (Alvarez et al. 2000), 
which results in the observed narrower distribution in I. We point 
out the overall similarity of the two AGB LFs, with only minor 
differences. 

The fainter peaks, on the other hand, show quite different be- 
haviour. The FWHM stays essentially constant with wavelength, 
with slightly narrower distributions for the SMC. Since the RGB lu- 
minosity function is increasing monotonically toward fainter mag- 
nitudes (Nikolaev & Weinberg 2000), the sharp cut-off in our data 
is a consequence of the falling detection efficiency in the com- 
bined OGLE-II+2MASS database. The dip between the brighter 
and fainter peaks is very close to the TRGB magnitudes determined 
by Cioni et al. (2000) from an analysis of about 150,000 objects in 
the DENIS catalog and that close proximity was interpreted both 
by Ita et al. (2002) and us in Paper I as the evidence for pulsating 
RGB stars. 

Following the referee's recommendations, we examined this 
agreement more closely. Due to the nature of the tip of the RGB, 
it manifests as an edge in the luminosity function (/obs)- Cioni et 
al. (2000) have shown very thoroughly in their Appendix that mak- 
ing an unbiased estimate of the discontinuity can be very difficult. 
Previous authors generally determined the position of the peak in 
fahs' while Cioni et al. (2000) argued for using f" bs . Adopting their 
considerations, we determined the second derivatives of smoothed 
luminosity functions in all bands (Fig. [3}- The OGLE-II LFs were 
smoothed with a Gaussian weight-function (FWHM=0.04 mag) 
and the derivatives were approximated by the difference ratios. 

The positions of maxima are listed in Table|2| For comparison, 
we also show the dips of the LFs and the TRGB values determined 
by Cioni et al. (2000). The overall agreement is excellent: most of 
the differences do not exceed the uncertainties. The largest discrep- 
ancy was found in Ks , which can, however, be explained by a ~0. 1 
mag systematic shift between 2MASS and DENIS ^-magnitudes. 
Therefore, we conclude the luminosity functions of OGLE-II vari- 
ables show exactly the same behaviour as, for instance, those of 
the much larger DENIS Catalogue towards the Magellanic Clouds, 
which included both variable and non-variable stars. 

Additionally, we find a significant wavelength dependence of 
the TRGB values, which follows exactly the changes observed in 
the tip of the RGB as a function of metallicity in globular clus- 
ters by Ferraro et al. (2000) and Ivanov & Borissova (2002). For 



instance, Ferraro et al. (2000) derived M^ p ~ — 0.6[Fe/H], which 
implies AK ~ 0.18 mag, assuming the metallicity differs by a fac- 
tor of 2 between the Clouds. The observed differences from the sec- 
ond derivatives are 0.19 mag, 0.25 mag and 0.23 mag in J, H and 
Ks, respectively (adopting A/io = 0.44 mag and correcting for 
differential extinctions) and they are also in good agreement with 
the relations in Fig. 4 of Ivanov & Borissova (2002). In / band, the 
difference is only 0.03 mag, which illustrates the insensitivity of 
AfxRGB (I) to metallicity and age (Lee et al. 1993). 

We interpret these agreements as our last piece of evidence 
that the dip between the AGB and RGB peaks indeed corresponds 
to the tip of the RGB. In Paper I, we have also considered the ar- 
guments of Alves et al. (1998) and Wood (2000) that all variables 
below the TRGB are thermally pulsing AGB stars and that their 
brighter cut-off matches the TRGB by coincidence. With the pre- 
sented behaviour, this explanation seems to be quite unlikely. Even 
if there are faint TP-AGB stars below the tip, their fraction must be 
small compared to first ascent red giants. 



3.3 Amplitude distributions 

We finish the comparative analysis by presenting the same six slices 
of the (period, amplitude, A's magnitude) data cube as were shown 
for the LMC in Fig. 4 of Paper I. Colours refer to three J — Ks 
colour ranges, enabling a rough classification of stars as "hot", 
"warm" and "cool" red giants. 

The distributions are similar to those of the LMC in that there 
is a good correlation between the amplitude and mode of pulsation. 
For amplitudes larger than 0.04 mag, practically all variables below 
the TRGB disappear. Furthermore, the colour distribution within 
any particular P-L ridge follows the same pattern as in the LMC. 
There are, however, some striking differences. Firstly, there is a 
relative lack of "warm" stars among the low-amplitude variables of 
the SMC. There is one blue dot for every 30 dots in the upper left 
panel; in the LMC, the ratio is 1:10. A similar under-abundance is 
also visible in the lower panels. Secondly, the lower three panels are 
dominated by the very red stars with J — Ks > 1.4 mag, which are 
probably carbon stars. For the highest amplitude variables (all with 
regular Mira-like light curves), there is only a handful of "hot" and 
"warm" giants, many of which fall below the fiducial Mira P-L by 
as much as 1.5-2 mag. These "faint Miras" be identified with the 
dust-enshrouded Mira stars (Wood 1998), since their extinctions 
increase strongly toward shorter wavelengths, reaching up to 4-6 
magnitudes in /. Our data illustrate nicely that, despite the fact that 
carbon-rich Mira stars obey the same P-L relation as the oxygen- 
rich Miras (Feast et al. 1989), the use of the Mira P-L relation in 
extragalactic distance measurements must be restricted to the bluest 
possible Mira stars. 



4 CONCLUSIONS 

In this paper we presented the results of a period analysis of OGLE- 
II red variable stars in the SMC. Multiple periodicity has been 
found for the majority of stars, and can now be identified as com- 
mon behaviour in pulsating red giants, regardless of the metallicity. 

We confirm the existence of distinct P-L relations below the tip 
of the Red Giant Branch, announced in Paper I for the LMC. Utiliz- 
ing the overall P-L distributions, we determined multi-wavelength 
relative distance moduli and extinctions, which support a consis- 
tent view that the differential distance modulus of the Clouds is 
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Figure 4. P-L relations in the SMC as function of the full amplitude of modes. Three J — K$ colour ranges were selected to plot in different colours. Turquoise 
(light gray): 0.9 < J — K§ < 1.2? blue (black): 1.2 < J — Ks < 1-4; red (dark gray): J ~ Ks > 1.4. This figure is available in colour in the on-line 
version of the journal on Synergy. Symbols in the lower panels are drawn larger for improved clarity. 



A/in m 0.44 mag. The luminosity functions revealed the AGB vari- 
ables follow very similar distributions in both Clouds. We showed 
that the dips between the AGB and supposed RGB peaks do in- 
deed correspond to the tip of the RGB (within a few hundredths 
of a mag). This is also supported by the wavelength-dependent 
AMtrgb (~0.20 mag in JHKs) found as an excess to the rel- 
ative distance modulus. It is consistent with the metallicity depen- 
dence seen among globular clusters and we take this as further ev- 
idence that the contamination by thermally pulsing AGB stars is 
likely to be small. 
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